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We show that the problem of stabilization of extra dimensions in Kaluza-Klein type cosmology 
may be solved in a theory of gravity involving high-order curvature invariants. The method suggested 
(employing a 3 slow-change approximation) can work with rather a general form of the gravitational 
action. As examples, we consider pure gravity with Lagrangians quadratic and cubic in the scalar 
curvature and some more complex ones in a simple Kaluza-Klein framework. After a transition to 
the 4D Einstein conformal frame, this results in effective scalar field theories with certain effective 
potentials, which in many cases possess positive minima providing stable small-size extra dimensions. 
Estimates made in the original (Jordan) conformal frame show that the problem of a small value 
of the cosmological constant in the present Universe is softened in this framework but is not solved 
completely. 
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Short title: Self-stabilization of extra dimensions 
PACS numbers: 04.50.+h; 98.80.-k; 98.80.Cq 

I. INTRODUCTION 



Modern cosmology is able to explain many properties 
of our Universe due to a substantial progress in experi- 
mental methods and technology and development of new 
theoretical models. Nevertheless, there is quite a num- 
ber of long-standing and challenging unsolved problems. 
We are going to discuss some of them. The first one, 
known for two decades, concerns the origin and form of 
the inflaton potential [jj] able to provide a sufficient pri- 
mordial inflation satisfying the observational constraints. 
Despite numerous attempts, see, e.g., the problem 

is yet to be solved. Another problem is related to the 
recently observed accelerated expansion of the Universe 
whose origin is usually ascribed to some unknown and 
invisible dark energy, which is likely to be identified with 
an extremely small but nonzero vacuum energy density, 
or cosmological constant. Certain problems of theoreti- 
cal origin are common to all modern theories, employing 
the idea of extra dimensions, among which well-known 
are M-theory, supergravities, diverse Kaluza-Klein and 
brane-world scenarios. A requirement which inevitably 
accompanies such investigations is that the extra dimen- 
sions should be invisible since our space-time looks 4- 
dimensional at least at scales over I0 -17 cm. In models 
where the extra dimensions are assumed to be compact, 
this is a severe restriction on their size. Moreover, to 
account for the observed constancy or nearly constancy 
of the fundamental physical constants (the gravitational 
constant G, the fine structure constant a and others 
4]), this size should be constant and stable, and there- 
fore a mechanism of stabilization of the extra dimensions 
is required in any scenario that involves them. 



The most popular way of reaching such a stabilization 
is to invoke a scalar field of non-geometric origin with 
an appropriate potential. This new entity seems some- 
what artificial because scalar fields already exist in this 
framework: the metric tensor components correspond- 
ing to the extra dimensions behave as scalar and vector 
fields in our 4D space. It is, however, difficult to obtain 
a suitable form of the potential for such fields 0, • To 
improve the situation, additional non-gravitational fields 
are invoked thus returning us to the starting point. 

On the other hand, the gravitational action, in gen- 
eral, includes terms containing nonlinear functions of the 
Ricci scalar as well as other curvature invariants. Such 
nonlinear terms inevitably appear if one takes into ac- 
count quantum phenomena Q, 0, Q • In fact, the first 
model of inflation was based on this idea [l(|. Grav- 
ity represents a wide choice of different forms of action 
since even renormalization of non-gravitational quantum 
fields against a curved background yields various curva- 
ture terms with indefinite coefficients, to say nothing of 
presently unknown classical "traces" of quantum gravity 
and Casimir-like contributions from a nontrivial space- 
time topology (if any). This wealth of possibilities has 
a negative imprint. Indeed, only one form of the gravi- 
tational action among many is realized in our Universe, 
and it is a challenging problem to choose the correct one. 

In this paper, we study curvature-nonlinear gravity in 
a space-time manifold with an arbitrary number of extra 
dimensions. Our aim is threefold. First, we intend to 
demonstrate that even the simplest version of nonlinear 
multidimensional gravity (with the Lagrangian F(R), R 
being the Ricci scalar) leads to a scalar field with po- 
tential and kinetic terms having promising shapes. Sec- 
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ond, we are going to point out some particular exam- 
ples of minima of the effective potential V with positive 
V = Vmin , which are, in principle, able to describe the 
present accelerated expansion of the Universe (certainly 
in a very rough approximation). Third, we try to find out 
whether or not, or to which extent, the well-known cos- 
mological constant problem can be solved in the present 
class of theories. 

In our study we use the slow-change approximation, as- 
suming smallncss of all derivatives, in a sense discussed 
below, and smallness of energy densities compared to the 
Planck scale. Such assumptions are rather widespread 
in studies of inflationary scenarios. This method allows 
dealing with rather general forms of the gravitational ac- 
tion and extended parameter ranges. With its help, new 
minima of the potential which can provide stability and 
smallness of the extra dimensions are found. 

The paper is organized as follows. In Sec. ^ we con- 
sider F(R) theories of gravity in multidimensional space- 
times where the extra dimensions comprise a constant- 
curvature space with a scale factor b(x) , depending on 
the external (observed) coordinates. Integrating out the 
extra dimensions in the action and applying the slow- 
change approximation, we obtain an effective scalar- 
tensor theory in 4 dimensions with a single scalar field 
b(x) and specific forms of its potential and kinetic terms. 
It should be noted that the complexity of the kinetic 
term can also affect the system dynamics and qualita- 
tively change its properties 0, . 

We also show that our approximation agrees with the 
conventional approach that reduces F(R) gravity to D- 
dimensional Einstein gravity coupled to a scalar field with 
a potential related to F(R) . 

Sec. IIIII is devoted to numerical estimates concerning 
possible stationary states of the extra dimensions. We 
briefly discuss the role of conformal frames, distinguish- 
ing the fundamental frame in which the theory is orig- 
inally formulated (saying what happens as a matter of 
fact) and the frame used to interpret the observations 
(showing what we can see) . The choice of the latter de- 
pends on the properties of references employed in our 
measurements ll2j . We here restrict ourselves to the sim- 
plest and maybe the most natural assumption, that the 
observational frame coincides with the fundamental (Jor- 
dan) one. Under this assumption, it turns out that the 
smallness of the present cosmological constant remains a 
problem though less severe than in conventional cosmol- 
ogy 

In Sec. lIVI we use our method for seeking stable states 
of the extra dimensions in quadratic multidimensional 
gravity, with F(R) = R + cR 2 — 2A where c and A are 
constants. This theory has been previously discussed in 
Ref. [l4|, where local minima of the potential V with 
Vmin < were found at values of R corresponding to 
dF/dR > 0. In the cosmological context, such minima 
can only lead to models with the anti-de Sitter space 
AdS 4 . We confirm this result but also find new minima 
of the potential in the unusual range where dF j dR < ; 



for some of them, Vmin > 0, which leads to de Sitter 
cosmology. 

Sec.|V]discusses other choices of the initial Lagrangian. 
Thus, with F(R) taken in the form of a cubic polynomial, 
we find minima of the effective potential with V m i n > in 
the region where dV j dR > . Then we demonstrate that 
our method can be successfully applied to a much wider 
class of multidimensional theories of gravity, e.g., those 
containing the Ricci tensor squared and the Kretschmann 
scalar. See l VII contains some concluding remarks. 

II. F(R) THEORY IN D DIMENSIONS 

A. Basic equations 

We consider a (D = do + d\ ) -dimensional manifold 
with the metric 

ds 2 = g^dafda? + e mx) b ab dx a dx b (1) 

where the extra-dimensional metric components b ab are 
independent of x^ , the observable space-time coordi- 
nates. The relevant case is certainly do = 4 but, in 
Eqs. we keep do arbitrary for generality. 

The D -dimensional Riemann tensor has the nonzero 
components 

RP ,v _ 

Jt pa ~ ±L pa i 

R»\a = B» := e-^e^), 

R- b cd = e-^R ab cd + S ab cd p^, (2) 

where capital Latin indices cover all D coordinates, the 
bar marks quantities obtained from g^ and b ab taken 
separately, /? M = d^(3 and 5 ab cd = 8 a c 8 b d - 5 a d 5 b c . The 
nonzero components of the Ricci tensor and the scalar 
curvature are 

Bp — + rf i B ^ 
R b a = e-^-R* a + 8 b a [UP + d x m% 
R = R[g}+ e~ 2l3 R[b} + 2d 1 nf3 + d 1 (d 1 + l)(df3) 2 , 

(3) 

where (df3) 2 = fi^ , □ = c^V^V* is the do- 
dimensional d'Alembert operator while R[g] and R[b] are 
the Ricci scalars corresponding to g^ v and b a t , respec- 
tively. 

Suppose now that 6 a j, describes a d\ -dimensional 
space of nonzero constant curvature, i.e., a sphere (k = 
1 ) or a compact d\ -dimensional hyperbolic space |l5j 
( k = — 1 ) with a fixed curvature radius ro normalized 
to the D -dimensional analogue m D of the Planck mass, 
i.e., ro = l/m D (we use the natural units, with the speed 
of light c and Planck's constant h equal to unity) . We 
have 

R a \ d = k m 2 D S ab cd , R b a = k ml (d x - l)5 b a , 

R[b] = kmldi(dx-l) = R b . (4) 
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The scale factor b(x) = e' 3 in JIJ is thus kept dimension- 
less; Rb has the meaning of a characteristic curvature 
scale of the extra dimensions. 

Consider, in the above geometry, a sufficiently general 
curvature-nonlinear theory of gravity with the action 



(5) 



S = \ m o~ 2 J V°(jd D x [F(R] + L 



where F(R) is an arbitrary smooth function, L m is a 
matter Lagrangian and g = | det (<7m iv ) | • The extra 
coordinates are easily integrated out, and the action is 
reduced to do = 4 dimensions: 



S= 1 -V[d 1 ]m 2 L 



gd 4 xe d ^[F(R)+L m ], (6) 



where 4 g = | det(g AU /)| and V[di] is the volume of a com- 
pact d\ -dimensional space of unit curvature. 



B. Slow-change approximation. The Einstein frame 

Eq. © describes a 4D theory which is nonlinear in cur- 
vature and, moreover, contains a non-minimal coupling 
between the effective scalar field f3 and the curvature. 
Let us simplify it in the following manner: 

(a) Express everything in terms of 4-dimcnsional vari- 
ables and (3(x); note that now 

R = Ri + <l> + h, 

Ri = R[g], h = 2diD/3 + d x (di + l)(9/3) 2 , (7) 

and we have introduced the effective scalar field 

0(a) = R b e - 2p{x) = kd, (di - l)m 2 e' mx) . (8) 

Recall that we have k = ±1 for positive and negative 
curvature in d\ extra dimensions, respectively, so that <f> 
has different signs in these cases by definition. 

(b) Suppose that all quantities are slowly varying, i.e., 
consider each derivative (including those in the defini- 
tion of R) as an expression containing a small parameter 
e\ neglect all quantities of orders higher than 0(e 2 ) (see 
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(c) Perform a conformal mapping leading to the Einstein 
conformal frame, where the 4-curvature appears to be 
minimally coupled to the scalar cf>. 

In the decomposition J7J), both terms /i and R 4 are 
regarded small in our approach, which actually means 
that all quantities, including the 4D curvature, are small 
compared to the D -dimensional Planck scale. So the 
only term which is not small is </> , and we can use a Taylor 
decomposition of the function F(R) = F(<fr + R4 + fi) : 

F(R) = F^+Ri+fx) ~ F(cf>)+F'( ( f>y(R 4 +f 1 )+..., (9) 



with F'(<fi) = dF/d(f). Substituting it into Eq. JBJl, we 
obtain up to 0(e 2 ) 

S = iv[di]m 2 J y/^d^x e dl/3 [F'(</>)i? 4 

+ F(</>)+F'(0)/i +L m ], (10) 

where (3 is related to (j> according to (JBJ) ■ The expression 
(|1(J|I is typical of a scalar-tensor theory (STT) of gravity 
in a Jordan frame. 

To find stationary points, it is helpful to pass on to the 
Einstein frame. After the conformal mapping 



6V ^ 9V = \f{4>)\9»v> 



F'W, (11) 



with the corresponding transformation of the scalar cur- 
vature 



lU = \I\{R i + jUf-^(dff 



(12) 



(the tilde marks quantities obtained from or with g^ v ) , 
the action (|10fl acquires the form 



S = -V[di] m 2 D I Vfjd 4 x{ [sign F'(cf>)} [R 4 + K] 



V{4>) 



-di/3 



F'(0) 2 ~ 

with the kinetic ( K ) and potential ( V ) terms 



(13) 



K = 3 (dff 



2 p 



2di- 



f 



(dpy 



v{4>) 



,-rfi/3 



F{4>)/F'{4>f 



(14) 



It remains to express everything in terms of a single 
scalar variable, say, <fi. We can write the action (|13(l in 
the form 



S 
L 



V[d 



1 2 



^(signF')L, 



Ri+ 2^Ein(0)W) 2 -y E in(0) 

L m = (signF')^ 7 ^L m ; 



(15) 



^Ein(</>) = TT-n 



1 

20 2 



iio "'(^7) -2di^ + ^i(rfi+2) 



Vkin(^) = -(signF') 



\0\m~L 



di(di - 1) 



F(4>) 
F'(4>f 



(16) 
(17) 



In I|12fl -lll5 [) . the indices are raised and lowered with g^ u ; 
everywhere F — F(4>) and F' — dF/d(f>. 

Eqs. lfTBTl - ((lTjl are valid for both positive (cf) > 0) and 
negative (<f> < 0) curvature of the extra dimensions. 
Minima of the potential VEin , defined in the Einstein 
frame, determine stationary points of the scalar field, 
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and they remain to be stationary points after a transi- 
tion to any other conformal frame, including the original 
Jordan frame. It is the Einstein frame that determines 
the scalar field behavior since its dynamics near an ex- 
tremum is governed directly by the potential V{4>) and 
only implicitly by the metric. 

If <j) resides at a minimum of the potential VEin j this 
potential turns into an effective cosmological constant. 
Being applied to the present cosmological epoch, it can 
determine the observable dark energy density that drives 
the accelerated expansion of the Universe. Alternatively, 
in principle, it may drive inflation in the early Universe. 

We have achieved one of our goals: the potential i|17|) 
valid for any function F{R) looks quite complex to have 
some nontrivial extrema. In addition, we have obtained 
rather a complex form of the kinetic term (|16[) . Its prop- 
erties are known to be as important for the field dynamics 
as the shape of the potential. For example, as shown in 
Ref. |13| . zeros and singular points of the kinetic term 
may be responsible for stable states of a scalar field. 

In our expressions, both the effective potential l|17|l 
and the kinetic term are singular at the values of <j) 
where /((/>) ~ F'(R) [the factor before R4 in the Jordan- 
frame action l|10|l ] is zero. Unlike many papers restricted 
to F'(R) > 0, we also include models with F'(R) < 0. 
As will be seen below, this opens new promising pos- 
sibilities such as new minima of the effective potential 
at which the extra dimensions may be stabilized. More- 
over, models with conformal continuations [19|,|20|, which 
unify regions with F' > and F' < , are possible; these 
regions correspond to different Einstein-frame manifolds 
but turn out to be smoothly connected in the Jordan 
frame. Though, such models are likely to be unstable, as 
follows from the experience of dealing with different so- 
lutions of scalar-tensor theories [2l|. Our main interest 
here is in other values of 0, namely, those at which it 
may be stabilized. 

In l|15|) - (|17|l we have actually changed the sign of the 
Lagrangian in the case F 1 < 0; to preserve the attrac- 
tive nature of gravity for ordinary matter, the matter 
Lagrangian density should appear with an unusual sign 
from the beginning. As a result, the sign of the whole 
action of gravity and matter will be unusual, without 
any effect on the matter equations of motion, and the 
conventional form of the effective Einstein equations at 
a stationary value of <p will also be preserved. Only the 
action of the <j) field itself can be unusual, according to 
Eqs. (|I6|I . HI7fl . It should be noted that the common sign 
of the total action does not affect quantum transitions 
as well. Indeed, the transition amplitude is expressed 
in the path integral technique as / exp (iS[q])Dq where 
q(t) is some dynamical variable. The transition proba- 
bility J J exp(iS[qi] — iS[q2\)DqiDq2 is invariant under 
the substitution i — > — i (with interchanging the integra- 
tion variables q\ and q^ ) . 



C. Comparison with the conventional approach 

There is a well-known conformal mapping [l6| bring- 
ing the theory (J5J to the form of general relativity with 
a minimally coupled scalar field \ that replaces the ad- 
ditional degrees of freedom connected with R : 



9ab = e- 2 "Wh AB , 



X = ± 



D-l 

D-2 



In \Fn 



J (X) = \F R \ 1/(D - 2 \(18) 



with Ffj = dF/dR. This results in a D-dimensional 
Einstein frame, with the action 



S = S DE [h] 



1 



n-2 



x \ (signF fl ) [R[h] + {d X )l] - 2V( X ) + L m e 
2V( X ) = e- D "(RF R -F), 



(19) 



where R[h] is the Ricci scalar obtained from Hab > h = 
I det(/i J 4s)| and (dx)1 = h AB x,AX,B', the initial Ricci 
scalar R is considered as a function of x- 

(An even more general transformation is known |17| . 
bringing a D -dimensional theory of gravity whose La- 
grangian contains an arbitrary function of two variables, 
f{R, <fi) , to an Einstein frame with two minimally coupled 
scalar fields, <fr and the one similar to x- Maeda's method 
|17| generalizes the corresponding mappings known for 
both scalar-tensor and f(R) theories of grav- 
ity) 

The theory (|19|l may be further reduced to 4 dimen- 
sions and brought with one more conformal mapping, 
now depending on (3(x) and similar to Ijlll) . to a 4- 
dimensional Einstein frame. Reduction in terms of the 
metric Hab gives 



S = S 4J [h] = l m ^ V [ dl ] J d^xV^h 



(signFfl) 



a := [3 + u, 



i? 4 N + (9x) 2 ,-di(di-l)(9a)2 



e- 2a R b 



2V{ X ) + L m e 



-Du. 



(20) 



where A h = \ det(/i M „)| . After a further transformation, 



(21) 



we arrive at a theory defined in the 4D Einstein frame 
with the metric 7 M „ . The action is 

S = S 4E [ 7 ] - im£- 2 V[d!] J ^d A x 

R[l] + (dx)l J + \d 1 {d l +2){da)l l 
- 2V( X ) e- dlQ + L m e- Du -^ a }, (22) 



< (signFij) 

lib 



_ e -{d 1 +2)a 



5 



with 7 = | detO^)] and {dx)* = J^X,nX,v 

Evidently, under the conditions of our slow-change ap- 
proximation, the action (1221) should lead to the same re- 
sults as 11131 . Let us confirm that it is indeed the case, 
assuming in (1221) that all derivatives <9 M involve a small 
parameter. We have, as before, © and similarly 



F R = F'{4>) + F"{<l>){R i + F l ) 



(23) 



where R± = R 4 [g] . Then, comparing the transitions 
9AB -* g^u -* g^u in dJ and g AB -> ->• /i/** -> 

7 M „ , we see that the two Einstein-frame metrics g and 7 
are related by 



Fr 



1 



F'{<t>) 



(Ri + h) 



7^ (24) 



and coincide in the leading order of magnitude. The lead- 
ing terms of the actions l|15fl and l|22|l . represented by 
their potential and matter terms, should then also coin- 
cide, and it is easy to confirm that it is really the case (the 
term RFr in (1221 transforms to cf>F'((j)) and is cancelled 
by the term with Rb , leading to the expression (|1TI) for 
the (j) ne ld potential). 

It is also seen that both fields x and a are reduced, in 
the same approximation, to the tf> field defined in Eq. JSJ: 



A, 



1 F" 
D-2F~>' 



D — IF 
D-2~Y 



(25) 

and — —0.^/(20). As a result, the whole kinetic 
term takes precisely the form <|16|) . So the action l|15|) is 
completely restored. 

One can note, however, that our approximation, being 
applied from the very beginning, not only makes things 
simpler but is more universal: it can deal with actions 
more general than (JSJ (see Sec.[Vj. 



III. EFFECTIVE CONSTANTS: SOME 
ESTIMATES 



Conformal frames and requirements to the 
model 



To relate the theory to Nature, it is necessary to spec- 
ify which conformal frame is to be confronted to obser- 
vations. 

In our problem setting, the Jordan frame described 
above appears to be fundamental since the initial theory 
is formulated in it, and the quantum effects that give rise 
to the nonlinearity of gravity are also supposed to take 
place there. 

It is, however, quite unnecessary to believe that this 
Jordan frame corresponds to the observed picture, in 
which the properties of our measurement instruments 
and the standards of the relevant physical units do not 
change in space and time The latter is thus the 

frame in which the basic atomic units (depending, above 



all, on fermion masses) are constant, if certainly they are 
all constant simultaneously. How sensitive is the matter 
Lagrangian to conformal mappings is even well seen from 
the coefficient before L m in Eq. I|15|) . to say nothing of 
a nontrivial metric dependence of matter fields entering 
into L m . 

Since we here do not specify how fermions are included 
into the theory, it is only possible to assume which is the 
observational frame. We will here dwell upon, probably, 
the most natural choice and make some estimates in the 
Jordan frame thus assuming that the basic and observa- 
tional frames coincide. 

As regards the Einstein frame in which the action takes 
the form l|13fl . we simply use it as a technical trick to 
determine stationary states. 

In any observational frame, in the expected station- 
ary state, the gravitational action has the approximate 
Einstein-Hubert form 



S, 



1 



GR 



gd 4 x{Rl -2A eff ), 



(26) 



where m 4 = (8irG N )~ 1/2 is the 4D Planck mass, Gn 
is the effective Newtonian constant and R\ is the ob- 
servable 4D curvature (for our choice, i?| = R4 ) . To be 
consistent with our approach and with observational con- 
straints, this stationary state should satisfy the following 
requirements: 

(i) A classical space-time description should be ad- 
missible, i.e., the true size of the extra dimen- 
sions should exceed the true D -dimensional Planck 
length 1 /m D , the fundamental length scale of the 
theory: 



m D b 



- /3o 



» 1, 



/m|<4(di-l). (27) 



(ii) The slow-change approximation should work, which, 
for (p = const , reduces to the requirement R4 <gC 



(iii) The observed size b* of the extra dimensions can be 
much larger than the Planck length I4 = I/TO4 w 
8 TCP 33 cm, but not larger than about 10~ 17 cm, 
which corresponds to the TeV energy scale. 

(iv) The predicted effective cosmological constant A e g 
should be very small to conform to the observa- 
tions: 



Aeff/m^ 



10 



-120 



(28) 



(the exponent is —120 instead of more conventional 
— 122 because the definition of 1TI4 involves 8w). 

Obtaining such a small value as H28|l is one of the well- 
known problems of theoretical cosmology since it is hard 
to explain without fine tuning why Km\ , generally asso- 
ciated with vacuum energy density, is so many orders of 
magnitude smaller than the characteristic energy densi- 
ties inherent to the known physical interactions (e.g., the 
Planck density m| for gravity). 
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B. Estimates in the Jordan frame 

Thus we assume that observations are performed in the 
Jordan frame. The corresponding action (|1(J|I is approx- 
imated by (|26[) with the effective constants 



8ttG 



JV 



V[di 



-2 id t pi 



A,-, 



F( 



2F> 



where 6 = e 13 ^ /m D = b 



(29) 

is the observable size of the 



extra dimensions and Fq — F'(4>o)- 

Eq. (|29() leads to the following relation between the di- 
mensionless quantities bom D and &0TO4 : 



i2 2 
On TO; 



V[d 1 ]^(6 TO D ) dl+2 . 



(30) 



The factor V[di] is of order unity; the same may be ex- 
pected from the dimensionlcss quantity Fq . By item (hi) 
,2^2 < 1Q 30 5 and Eq ^ jgg gives 



above, 65m 



ooTO D < 10 30 /( dl+2 ), 

m 4 /m D ~ {b m D ) dl < i0 15( V( d i+ 2 ). ^ 

Thus Tti D not too much differs from vri\ . It means, in 
particular, that our slow-change approximation (item (ii) 
above) works manifestly well in almost all thinkable cir- 
cumstances since R4 is the observable curvature. 

If m^/riio is sufficiently close to unity this approxima- 
tion is even valid for the curvature characteristic of pri- 
mordial inflation at the Grand Unification scale: thus, an 
a priori estimate is R^/m 2 ^ ~ 10~ 6 . In particular models 
this inequality is strengthened. 

Thus, at chaotic inflation governed by the inflaton tp, 
the Ricci scalar is expressed in terms of the Hubble pa- 
rameter H and the potential V(ip) as 



R A = 12H 2 : H = - 



inG 



(32) 



where the term with ip is omitted since if is assumed to 
be slowly rolling down the slope of the potential. Let us 
estimate R4 in the quadrativ model of chaotic inflation, 
The observational data indicate that 



with V = \m 2 (p 2 



H: 



1(T 6 TO4 



during inflation, and so 



R A /m\ = l2H 2 /m 



10" 



(33) 



(34) 



Due to the slow rolling, one may not worry about the 
smallness of fx at the inflationary stage. The kinetic 
terms are, however, the largest at the end of inflation, 
at reheating, when ip ~ 7714 , while the slowly changing 
Hubble parameter becomes comparable with the inflaton 
mass, H ~ m, i.e., m/m^ ~ 10~ 6 . Reheating is a result 
of quick oscillations of the inflaton around the minimum 
of V , where ip ~ rrvp , hence here we also obtain a small 
parameter: 



ip/m] 



10" 



(35) 



which justifies our approximation provided TO4/to d < 
10 6 . Otherwise our approximation starts to hold at later 
stages of the evolution, and it is manifestly good for the 
present-day Universe. 

Let us now return to the inequalities l|31() and estimate 
A e ff. The ratio A e g/m 2 may be expressed as follows: 



A 



cir 



F( 



,2 u2 



2m 2 D F> 
1 F(<f> 



-2d! 



2*2 



■{F^ ) V[d 1 }) d i+^ (TO 4 6 ) d i+ 2 . 



(36) 



If we try to make this ratio small, the last factor in 136|) 
can give at best about 10 -30 (for sufficiently large d\). 
The remaining 90 orders of magnitude must be gained 
due to unnatural smallness of the dimensionless quantity 
F(4>o)/m 2 - ) (actually of the initial cosmological constant) 
and/or greatness of Fq. The latter variant is, however, 
unsuitable since it will make impossible Oq.to d ^> 1 , see 

As a result, the problem of fine tuning remains topi- 
cal, though it should be noted that the very small value 
of 10 -30 appears in this approach without any artificial 
effort. 



IV. QUADRATIC GRAVITY WITH A 
COSMOLOGICAL CONSTANT 

In what follows we consider pure gravity ( L m = ) and 
use the units m D = 1 , thus dealing with dimensionless 
quantities. 

As the first and simplest example of multidimensional 
nonlinear gravity, consider (|SJ) with the function 

F(4>) = (t) + C (j) 2 - 2A, c, A = const. (37) 

Then Eqs. 1|16|) and l|17fl give the effective potential 

j-di/2 



v£iW =-sign(l + 2af>)[d 1 (d 1 - 1)]" 
/2 c0 2 + ^2A 
m (\ + 2c<t>) 2 

and the coefficient ATEin of the kinetic term 



(38) 



= 



1 



(1 



2c(t>) 2 (j> 2 

j2 



c 2 (j> 2 {d\ - 2dx + 12) 



dfcq 



2) 



(39) 



In Ref. the authors studied the case F'(R) > 
and showed analytically that a minimum of the po- 
tential stabilizing the extra dimensions is only possible 
with a negative energy density and a negative curvature 
k = — 1 . We confirm this result in our approach, see 
examples in Fig. 1. The energy density corresponding to 
the minimum value of the potential is negative, which, in 



FIG. 2: The same as in Fig. 1 but with c = 1.5. The minima 
of V(<f>) are situated in the range <f> < —1/3, where F' < 0. 



FIG. 1: The potential Vbin(</>) for di = 3, c = 0.3 in the 
range <f> > —5/3, where F'(R) > 0, for a few values of A 
(indicated near the curves) . The dashed line ( A = 1 ) shows 
a minimum at = 0. The values of if are expressed in the 
units m D = 1 . 



the cosmological context, corresponds to anti-dc Sitter 
space-time. 

There is one more minimum of the potential Vein 
with F' > 0, existing in the range A > and located 
at the point = 0. The asymptotic 0^0 corre- 
sponds to growing rather than stabilized extra dimen- 



sions: b = e^ ~ 1/y |0| — > oo . A model with such an 
asymptotic growth at late times may still be of interest 
if the growth is sufficiently slow and the size b does not 
reach detectable values by now. 

Let us check whether it is possible to describe the mod- 
ern state of the Universe by an asymptotic form of the 
solution for — > as a spatially flat cosmology with the 
4D Einstein-frame metric 



dsl = dt 2 - a E 2 (t)dx 2 



(40) 



where o E is the Einstein-frame scale factor. Two inde- 
pendent components of the Einstein-scalar equations for 
(3(t) and a E (t) are 



(3 + 3^(3 



2A 



a E di + 1 
3^ = U(rfi + l)/3 2 + A e-W. (41) 



We seek their solution at large t in the form 

a E (t) oc t p , e 13 w (t/t*) q , p, q, t* = const > 0, 



and find that such a solution docs exist with 



(42) 



p = 3h, q = 2/d\, 

h:= + 



Passing over to the Jordan frame with 

ds\ = dr 2 - a 2 (r)dx 2 = jdsi, (44) 

due to F'(0) w 1 , we can put simply / = c dl13 w (i/i*) 2 . 
We obtain 



t/U oc c T ' u , o(t) oc e ^h-i)r/u t fe(r) K e : 



Ktl = h 2 {9h- 1). (43) 



2r/(dit») 

(45) 

The external scale factor o(r) grows exponentially, 
which conforms to modern observations if one properly 
chooses the constants. The internal scale factor 6(r) 
grows much slower for sufficiently large d\ , but the vol- 
ume factor ~ b dl grows approximately at the same rate 
as a(r), which means that, e.g., the effective gravita- 
tional constant will change too rapidly, at nearly a Hub- 
ble rate, contrary to observations. We conclude that the 
model with — > at late times is not very promising. 

We have also tried to find minima of the potential in 
another region, where F' < 0, in a wide range of the 
parameters c and A. 

Some examples of the behavior of the potential for the 
specific value c = 1.5 are shown in Fig. 2. Potentials with 
A < show minima with V m i n > 0, so that A e a > 0, 
which should lead to cosmological models with de Sitter 
external space and stable extra dimensions. Such models 
again correspond to < 0, i.e., hyperbolic extra dimen- 
sions. By fine tuning of the initial constants c and A 
we can obtain the present-day values of the cosmological 
parameters. 

As in Fig. 1, we also find a minimum at = 0, whose 
properties have already been discussed. 

We see that the behavior of the system is drastically 
different in different ranges of and depends on the nu- 
merical values of the initial parameters. We also confirm 
that gravity alone can stabilize the size of extra dimen- 
sions, without need for introducing other fields with spe- 
cific forms of potentials. Positive values of the effective 
cosmological constant A e fj (vacuum energy density) were 
found in the range where F' < 0. This is a feature of 
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FIG. 3: The potential Vkin for cubic gravity with d\ — 3 and 
c — 2. The curves are labelled by values of the parameter C . 
Note that F'(<f>) > near the minima. 
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FIG. 4: VBi n (<?i) for cubic gravity with di = 3, c = — 1 . Lines 
are labelled by values of the parameter C . Again F'(4>) > 
near the minima. 



quadratic gravity only: as will be seen in the next sec- 
tion, e.g., for cubic gravity one may have A e s > where 
F' > 0. 



OTHER THEORIES 



A. Cubic gravity 

The cubic theory of gravity with 

F(R) =R + cR 2 + CR 3 



(46) 



is another example of curvature-nonlinear gravity for 
which Eqs. i|16|) and (|17fl may be applied. The shape of 
the potential appears to be very sensitive to the values 
of the parameters c and C . This property may be used 
to adjust the potentials in such a way that its minimum 
in the Einstein frame will lead to an extremely small but 
positive value of A e g- in the Jordan frame, and it turns 
out that, unlike quadratic gravity, this can be achieved in 
the range where F' > . Examples of such a behavior are 
shown in Fig. 3. Thus cubic gravity can in principle lead 
to an appropriate description of our Universe with com- 
pact and stable extra dimensions. At large times, when 
the influence of matter may be neglected compared to the 
effective cosmological constant, we shall obtain de Sitter 
expansion of the three observable dimensions. 

The range of positive curvature, <f> > , appears to 
be of lesser interest. It also contains minima but only 
with Vmin < 0, see Fig. 4. The states in such minima 
are metastable (since the minima are only local), so that 
such a space would have a finite lifetime. The latter, 
however, could be made arbitrarily long by choosing the 
parameters. 



B. Extensions 

The method discussed above allows considering wider 
classes of Lagrangians. Let us briefly demonstrate 
this by adding terms proportional to the Ricci tensor 
squared RabR AB and the Kretschmann scalar JC = 
RabcdR ABCD ■ By common views, these and other 
high-order curvature terms appear due to quantum cor- 
rections, and it seems natural to include them on equal 
footing with cR 2 . Now the action has the form 

S = ¥M f d 4 x y/^e d ^(R + cR 2 
2k 2 J 

+ Cl R AB R AB + c 2 K - 2A) , (47) 

where the internal variables have been integrated out in 
full analogy with Eq. JSJ) . For the metric , it is easy to 
obtain expressions for R AB R AB and K.: 



RabR ab = R^vR^ 1 



IdxR^B*- 



i K a bK , 



+ 2 e- 2p R[h] [130 + di (a/3) 2 ] + di [□ p + dx (dp) 2 } 2 

(48) 

K, = JC[g] + AdxB^B^ + e" 4/3 /C[/i] 

+ 4e~ 2/3 i?[/i](d/?) 2 + 2di(di - l)[(d/?) 2 ] 2 . (49) 

In the slow-change approximation, in the same manner 
as in Sec.[n] we obtain the 4D effective Lagrangian 



gL 



Ri(l + 2ccj>) +</>+ c tot cf> 2 - 2A 

(l + 2c0)/i +2ci0D/3 + 2(c 1 d 1 +2c 2 )(9/3) 2 

(50) 



with 



Ctot = c + — 



2c 2 



g?i d\(d\ — 1) 
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The conformal mapping leads, after some calcula- 
tions, to the Einstein-frame Lagrangian i|15[) with the 
kinetic and potential terms 



c\ + c 2 
2<j)(l + 2ccj)) : 



V E in(4>) = -sign(l + 2c0) 



- 1) 



di/2 



ctot4> 2 + (t> ~ 2 A 
(1 + 2c0) 2 



(51) 



> (52) 



where the term K^ n ((f>) is taken from (|3l?l) . 

The presence of the parameters c\ and C2 adds free- 
dom in choosing the shape of the potential. The kinetic 
term also acquires a more complex form which could sig- 
nificantly affect the field dynamics. Thus, as shown in 
Ref. |13| , zeros of the kinetic term can represent station- 
ary values of 4> ■ R means that a field could be captured in 
the vicinity of such points in addition to minima of the 
potential. An analysis of kinetic terms like (|51|l could 
lead to possibilities of interest, and we hope to return to 
this point in our future work. 



VI. CONCLUDING REMARKS 

We have shown that multidimensional gravity alone, 
even without any fields of non-geometric origin and with 
a very simple choice of the geometry, can be a basis for 
rather complicated phenomena. The assumption on non- 
linearity of gravity, used in this study, inevitably follows 
from quantum corrections to the Einstein gravity and 
must be taken into account for completeness. 

We have introduced the slow-change approximation in 
nonlinear multidimensional gravity, making the analy- 
sis easier. This approximation proves to be valid for a 
wide variety of phenomena where the curvature and en- 
ergy scales are far from Planckian (though Planckian in a 



multidimensional sense, which may mean scales different 
from those known in four-dimensional physics). 

Using some simple examples of nonlinear gravity with 
an arbitrary number of extra dimensions, we have ob- 
tained an effective scalar field with quite a complex form 
of the potential. The kinetic term is also nontrivial and 
adds complexity to the effective field dynamics. The po- 
tentials possess minima, both stable and mctastable, thus 
stabilizing the size of extra dimensions. Some of them 
seem quite promising for the description of the present 
state of the Universe with a small positive effective cos- 
mological constant and stable and small enough extra 
dimensions. 

A problem that could not be solved in this paper is that 
of choosing the observational conformal frame, which is a 
necessary step in confronting the theory to observations. 
For our numerical estimates we identified the initial (fun- 
damental) and observational frames. One should, how- 
ever, bear in mind that it is not the only possible choice, 
and the right one should follow from a full underlying 
theory. 

The precise number of extra dimensions d\ was not 
very important in our study, and we used the only value 
d\ = 3 in numerical calculations, though the method 
is valid for arbitrary d\ . This number may prove to be 
more essential in more complex geometries and nonlinear 
gravity theories. 

One of possible applications of this method is the brane 
world scenario, where, in particular, it opens a direct 
way of obtaining a stabilizing radion potential instead of 
postulating it. 
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